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Analysis of Simultaneous Skylab and 
Ground Based Flare Observations 


1 . INTRODUCTION 

The purpose of this program has been to reduce and analyze Eel and 
Hell resonance line data from Skylab and make comparison mth Hel line 
intensities taken simultaneously from the Lockheed Eye Canyon Solar 
Observatory. Specifically, ve will obtain relative and absolute intensities 
where possible as functions of position and time ^rithin the most interesting 
solar flare observed simultaneously from Skylab and the ground. The observed 
intensities will be compared with theoretical calculations made as part of 
this study with the objective of inferring electron temperatures and 
densities within the flare as a function of time. 

' The investigation described here has been a collaborative effort 
between us and the Skylab NRL principal investigating group. Ne would like 
to express our appreciation to Ik*. R. TOusey for providing data before pub- 
lication. Dr. J- D. Bohlin provided us with flare lists and much advice and 
Dr. Fred Rosenberg generously supplied tapes with microdensitometer informa- 
tion on many plates taken with the NRL XUV Spectroheliograph. This data is 
discussed in Sec. II. C. 

Hel D^ filtergrams of the solar disc were obtained from May 1973 to 
February 197^ at the Lockheed Eye Canyon -Solar Observatory. These filtergrams 
contain the emission Tr/ithin an 0.4A. band about the center of the line 
( 5376 A.) and were taken at 15 or 50 second intervals for most of this period. 
The purpose of such patrol observations was to catch active 
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regions in early developing stages. The filtergram field of view is 8’ x 10’ 
while the spatial resolution is of the order of 1”. The filtergrams thus 
contain a great deal more spatial ixiformation than corresponding spectrographs. 
The data is discussed in Sec. II. B. 

Recent data show many complex and as yet unexplained phenomena in the 
line (Ramsey, 1970; Harvey, 1971; Ramsey, et al. , 1975)* 

The emission has shown intself to have the potential of providing signifi- 
cant new information on dense flare centers where its emission is generally 
concentrated. The following general conclusions have heen drawn concerning 
flare structure. is observed in emission and/or absorption against the 
disc during most flares of importance 1 or greater. emission is generally 

confined to larger disc flares and occur? more frequently in limb flares. The 

\ 

emission component of D, flares corresponds to the brightest and spectrally 
broadest parts of the H - a flsire. The localized bright centers in are 
much shorter lived and usually show more rapid changes than can be seen in 
H - a filtergrams. The brighter emission centers are visible for tens of 
minutes. In most cases the e m ission elements were adjacent to or over- 
lying sunspots. A wide variety in flare structure has been observed. 

absorption elements correspond generally to weaker parts of the 
3 

H - a flare and are longer lasting than emission. The absorption components 
of flares aside from filaments, surges and loops, are not visible in the 

red wing at + 0.8 a. About half the emission components sho^r an initial 
red wing brightening usually lasting for .about 30 sec. This is interpreted 
in light of various observational evidence as an initial downward motion of the 
emitting region. 
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For purposes of comparison with Lockheed D, data it was decided that 

3 

the I®RL experiments could provide the most complete and useful complimentary 
information. The HRL instruments on Skylab were the extreme ultraviolet 
spectroheliograph (so 82A) and the ultraviolet high-resolution spectrograph 
(S082B) . The S082A instrument is a slitless spectrograph having a short 
wavelength range 150-535A and a longer range 32I-65OA. The images of the He 
resonance lines are fairly well separated, giving total Integrated line 
intensities. S-O82A images the entire solar disc for all wavelengths in rapid 
time sequence. The spatial resolution for the He resonance lines is estimated 
at 2". Time resolution depends on the exposure time, which varies from the 
shortest exposure of 2.5 sec. to a maximum manual exposure of 48 min. The 


SO82B instrument also has two wavelength ranges (970-19701 and 1940-59^1) • 
Spectral resolution is .04 - .08I with a field of view of 2" x 60". There is 
no spatial resolution along the slit. 


We have developed computer codes under the SOAP program which allow 
for the calculation of total He line intensities and line profiles from model 
flare regions. A complete description and results are given in our final 
summary report (Kulander, 1973) which we shall refer to as Rl. These codes 
incorporate simultaneous solution of the line and continuum transport equations 
as needed together with the statistical equilibrium equations for a 3O level 
Hel, Hell, Helll system. The geometrical model is a plane -parallel layer 
irradiated on one side hy the photospheric radiation field. A statistically 
steady state and uniform electron temperature and density with position were 
assumed. The energy level model consists of all terms through principal 
quantum number four.. Our study has been confined to conditions we believe 
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characteristic of He flare regions, namely electron temperatures between 10 

ijQ PQ pij. 

and 5 X 10 K and electron densities between 10 and 10 cm Our codes 
also allow the solution of the statistical equilibriiun equations for arbitrary 
• values of the various line and continuum radiation fields. Results of the 
inteipretation of the observed data in terms of these parametric solutions and 
with simultaneous solution of the transport equations are discussed in Sec. III. 
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II. A. Selection of 


f iltergrams were taken at Lockheed during the period May 28, 1973 to Feb. 4, 1974* 
Generally the regions studied coincided with active regions under observation by 
Skylab. Photographs were taken every 15 or 30 sec. depending on seeing conditions. 
This data with observation times, seeing condition^ etc. , is sujnmarized in the 
World Data Center A Report UAG-43, Coffee^ 1975# 

A list of flares observed by the NRL S082A instrument was compiled by J. D. Bohlin 
and V. E, Scherrer and generously supplied to us by Dr. Bohlin. This list is given 
in Table II. A. 1. A total of 17 flares was observed, of which about a dozen had 
reasonably complete coverage. The list is a result of the merging of flare logs 
compiled independently from two different approaches. One log was made from a 
plate -by -plate inspection of the entire mission, using as a criterion whether or 
not He II 304 was solarized and the presence of emission lines from the high stages 
of ionization (primarily in the short wavelength band) , The other log was com- 
piled by comparing 508 2A coverage with the SOLRAD data. 

Of these 17 flare events we have simultaneous coverage for only 7. These are 
the flares of 15 June 73, 7 Aug. 73, 9 Aug. 73, 31 Aug# 73, 1 Sept. 73, 5 Sept. 

73, and 7 Sept. 73. Our films for these seven flares were studied in detail. 

For two of these events only two frames were taken by S082A. Since as we shall see 
shortly several exposures at each time are normally required for good intensities, 
these events would not yield much time dependent information. The 15 June 73 event 
occurred very early in the morning in Los Angeles and our data is of very poor 
quality. The 1 Sept, event is very close to the limb and the contrast is poor 
on our film. The 5 Sept, event shows minor activity 'in and our coverage does not 
begin until after the maximum. For the 7 Sept, event neither NRL nor Lockneed has 
coverage at flare maximum. 

Thus’ the most promising event where there 'is both S082A and D^ coverage through the 
flare maximum is that of 9 Aug. 73. Flare xaaxiinum is near 15:53 UT. NRL has 14 
frames between 15:20 and 15:58 while we have D^ filtergrams every 30 secs, under 
good seeing conditions. 
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TABLE II. A. 1 
NKL S082A FLARE LIST 


Date 

NBL 

Plate 

Nos. 

Time 

(NT) 

1 Position 
i Angle 
(+5°) 

Notes 

15 June 1973 

1A312- 
I 347 

14:11- 

14:38 

N20W45 

1 

WV LG & ST. M2: Max. 14:14. 
Complete sequence. 

7 Aug- 1973 

2A006- 

007 

18:47 

18:49 

N10W30 

C3: Max. 18:47- Both plates LG W. 
Many flare images superposed on film 
defect. Exposure one hour earlier and 
later still show AR to be very bright. 

9 Aug. 1973 

2A022- 

035 

15:20- 

15:58 

N10W50 

M2: Max* 15:53* May be sequence of two 
flares, with small limb flare on 2A022 
the cause of large EPL. Many flare 
images superposed on film holder defect* 

31 Aug, 1973 

2A356- 

357 

21:57 

21:59 

1 

S05W15 

Two exp. (WV ST.) of small flare. 
Emission lines well- separated and 
compact . 

1 Sept. 1973 

2A371- 

373 

21:26- 

21:32 

N05E80 

Subflare. Emission lines not distinct. 
LG 6c ST WV. Max. 21:15. 

2 Sept. 1973 

2A376 

00:46 

S10E70 

Subflare. Emission lines not distinct* 
ST WV. Max. 00:43 

4 Sept, 1973 

2A398- 

399 

m 

S20E70 

Txiyin solarized spots in He II. Bright 
footpoints in many highly ionized lines 
OT LG 6c ST j synoptic pair. Max. 11:30 

4 Sept. 1973 

2A403 

16:36 

NIOEIO 

C6: Max. 16:35 

5 Sept. 1973 I 

2A419- 

420 

18:31- 

18:37 

Disk 
Center • 

Distince arch visible 
ST WV only. 

5 Sept. 1973 

2A42X- 

422 

■1 

S20W30 

Max* @ 20:10. Small’ flare in NE quadrant. 
Complex emission structure 

7 Sept. 1973 

! 

2A432- 

448 

12:21- 

18:28 

S20W50 

! 

1 

1 

XI flare w/4 consecutive peaks. Post flare 
loops over neutral lines esp. apparent. 

ALl ST OT. Max. 12:10. 
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Date 

NRL 

Plate 

Nos. 

Time 

(TJT) 

10 Sept. 1973 

2A459- 
; 460 

02:32 

2 Dec. 1973 

3A037- 

045 

15:04- 
i 15:17 

16-17 Dec. 73 

3A097- 

143 

19:00- 

02:00 

22 Dec. 1973 

3A176- 

178 

00:22 

01:20 

15 Jan. 1974 

3A375- 

14: 25- 


376 

14:28 

21-22 Jan. 74 

3A452- 

472 

23:15- 

00:49 



S20W60 


S10W80 

S20E85 

S10E05 

N10W85 

N30W50 


Notes 


Synoptic pair may have caught small 
flare. He II barely solarized. LG 
& ST WV. 

Ml flare. Max. 15:18 LG & ST WV. 

Good sequence. 

Complete sequence on C2 & Ml (Max. 
00:32 & 00:41) limb flares. 47 plates 
total. 

Small flare on 3A176 (ST WV); emission 
lines still apparent on synoptic pair 
177 & 178, May only be a very bright 
AE. 

Limb flare (C6: Max. 14:25) with 
sheet- like surge, ST WV only. Surge 
seen hour earlier as well, 

C8 flare, full sequence of 21 plates. 
Max. 23:23, 
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IX, B. Filter gram Data Analysis 

The flare event of 9 August 73 occurred in the active region designated 185 
quadrant) on the NOAA daily charts. Fig. Il.B.l shoxjs the NOAA schematic of the Sun 
at 14:20 UT. The region was observed in the line at Rye Canyon for several 
days before and after the flare. As originally observed the spot group was very 
complex with many small and intermediate size spots. By 9 August only the largest 
spot and one of the smaller eastern spots remained. The spot centers are about 
1-1/2’ apart and are located roughly at position angle N 10^, W 50^. The smaller 
Spot lies roughly E of the larger one. The central intensity of the Larger spot 
Is as low as 357o of the continuum while in the smaller spot it is about 75%, 

Figs, XI, B, 2, 3, and 4 are photographs of the spot region through the filter 
at 15:46:00, 15:53:30 and 15:57:00. Exposure times were a small fraction of 
a sec. The first photograph shows the pre-flare undisturbed environment, the 
second photograph the maximum and the third photograph the decay phase, The W 
limb is visible in each exposure. 

The flare or flares occur in the smaller easterly spot. Visual inspection of the 
film shows two distinct bright emission points on each side of this spot. These 
two bright flare kernels can best be seen in Fig, XI. B. 3 at 

15:53:30, Both kernels lie roughly in a lil - S line. The northern kernel is first 
visible at 15:52 fading out at 15:58 while the southern kernel is first visible 
at 15:53 and fades out at about 15:57, Both kernels appear very bright against 
the spot background. They appear to represent two distinct events at different 
locations , 

Two sequences of densitometer tracings have been made from our original 8 -mm 
film. The first consists of tracings of the four best coincident frames 
listed, in liable II. C, 2 at 15:46:00, 15:53:30, 15:54:00, and 15:57:00, The 
slit size used was 8x8 microns. Digital output for intensity was obtained using 
our H&D curves for 80 x 80 micron and 32 x 32 micron averages. 32 microns corre- 
sponds to 1,0” arc. The area traced was 106” x 160” and is shown in Fig. Il.B.l, 
Figs, IX. B. 5 and 6 show Intensity contours of the smaller spot from the tracings 
at 15:46:00 and 15:53:30, The two distinct emission areas are prominent at 
15:*53:30, The area covered by the Figs, is 25” x 30”, The contours are labeled 
with the percentage of the continuum intensity. 
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Figure II. B,:? 


Dj Filtergram Intensities in Percent of Nonnal Continuum. 

TJI 15:^16:00 
Area 25" x 29" 
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second sequence of tracings represents a similar area for each of the 
tvrelve frames taJ£en every 30 sec* beuvreen 15:51*00 and 15 "56^^50 tJT* The slit 
size used at this time "was 24 x 24 microns or 0*75”« Intensity contours for 
these frames are shoim in K-gs. II*B*7 to IIoB*l8, The time and average back- 
ground intensity are given on each figure. The contours represent 3^0; 400^ 

440^ 480 and 520 arbitrary intensity units. Each figure shows an area of 
approximately 23" x 34". The contours are labeled mth the relative intensity 
values. The N and W direction is indicated on each figure. Areas Trdthin the 
spot with intensities greater than the normal continuum, are shown on a number 
of the frames. There appear to be many instances of very small bright emission 
points that last for only one frame; see for instance Figure H.B.8. The 
southern flare kernel is first seen at 15:52:00 mth the northern flare kernel 
beginning at 15:52:50* This is followed at 15 "5 3 ^^00 by a slight diminishing in 
intensity in both flares (possibly due to the initial red shifting noted by 
Eamsey; 1975) followed by the greatest intensities at 15:53:30 in the northern 
flare. 

The northern flare becomes much more intense than the sour them flare. The maximum 
intensity at any one point relative to the undisturbed continuum is 1.26 in the 
northern kernel and 1,10 in the southern kernel. With respect to the underlying 
spot backgrounds these maximum relative intensities become 1,62 and 1,42, The 
area of the northern region becomes much larger at maximum while that of the 
southern region remains about the same. The intensities in both regions slowly 
decay still being slightly above background at 15:56:30. These features are 
illustrated in Figs. II, B. 19 and 20 where we show the intensities at four dif- 
ferent locations within the flare for each frame relative to the initial spot 
intensity at that location. 

We note the rapid rise to maximum followed by the slower decay. The southern kernel 
reaches maximum intensity at 15:52:00. *In bothflares there seems to be an oscilla- 
tion in intensities with a period of 1 - 1-1/2 minutes. 
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Figure H.B.17 Relative D_ Filtergram Intensities 









Figure H.B.19 


Intensities Relative to Underlying Spot 
for h Positions in South Flare Kernel. 


Background 
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II • C. SO82 Data Analysis 


A list of the 14 URL plates for the flare of 9 Aug. 73 is given in Table II.C.l 
along with exposure times, wavelength range covered, and comments. We can see 
that only the four plates 23, 29, 30, and 31 encompass the wavelength range for 
the resonance lines of Hel, i.e., the long wavelength range. Plate 23 is vastly 
overexposed while we note that on plates 29, 30, and 31 the 584 X line is so 
badly out of focus as to render quantitative analysis impossible. The longer ex- 
posures on the short wavelength plates give a solarized flare kernel at 304 S 
limiting analysis to exposures of 10 sec. or less. Plate 32 is fogged due to pres- 
sure of the backing plate against the film. The usable plates for our purposes 
are then 24 - 27, 33, and 34. These plates will only provide information on the 
resonance lines of Hell. 

Since plates 26 and 27 and plates 33 and 34 are taken at essentially the same time 
we chose the shorter exposure time plates 26 and 34 for analysis. The four most 
useful frames ^ together with the closest filter gram times are given in Table 

II. C. 2. There should be no large error due to different exposure times. 

r* 

Fred Eoseriberg at IJRL has supplied us with tapes containing digitized output 
of mierodensitometer scans of S082A data. Plates 20, 2k - 2&, 32, 33 and 55 
were scanned for the Hell 30k, and 2k3k images . Itie H & D CTxrve was 
provided for the 504 a line and it was used also for the 256 and 245A lines. 

Some of these plates were considered for their 256 and 2^\. intensities even 
though the 30^A results are solarized. Most of the scans were made of a 
2500 X 2500 micron region as shown by the dashed lines in Fig. II. A. 1. The 
slit size used was 10 p. x 10 p,. With a solar diameter of I8.5 mm. the 10 p 
steps correspond to a 1.04" resolution. 

We have unpacked the tapes sent to us by MEL and have reconstructed the output 
and calculated intensities for 20 p or 2.08" average slit size. 304 a total 
line intensities were obtained for plates 2 k, 23 , 26 , 32 and 54. 256I and 

243a total line intensities were obtained from plates 24, 26 , 28 and 52. Pigs. 
II.C.l - II. C. 3 show the sun at 304 a from plates 24, 26 and 54. The flare 
is prominent at the eastern edge of AR 185 in all three photographs. 
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TABLE ll.C.l 


SO 82 Plates For 9 August 1975 Flare 


Time 

TJT 

Plate 

No. 

Wave 

Length 

Exposure 

Time 

(sec) 

Comments 

15:20:32 

2A022 

ST 

40 

OK but exp is long for Hell. EPL 
in progress j possibly from small 
limb flare. 

15:35 

023 

LG 

160 

Wrong roll, and exp. too long 

15:46:07 

024 

ST 

9 

Looks good. Exp. short enough so 
Hell is not solarized 

15:30:20 

025 

ST 

1 3/4 

OK. exp. very short; flare kernel well- 
defined 

15:54:00 

026 

ST 

2% 

OK. exp. very short; flare kernel well- 
defined 

15:54:05 " 

' 027 

ST 

\ 

\ 

10 

OK.*flare kernel is solarized but still 
well-defined 

15:54:17 

028 

ST 

40 

Exp. much too long for Hell -solarized 

15:54:58 

15:55:05 

029 

030 


LG 


He I 584 out of focus. 

15:55:27 

031 ' 





15:56:49 

032 

ST 

2% 

Flare fogged by plate defect 

15:56:53 

033 

ST 

10 

Flare slightly solarized 

15:57:05 

034 

ST 

% 

Flare exposed OK 

15:58:20 

035 

ST 

40 

Flare solarized 


ST =: Short Wavelength 
IjG- = Long Wavelength 
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TABLE II. C. 2 


COINCIDENT FRAME TIMES 
FLARE OF 9 AUGUST 1973 


NRL 

Plate 

Number 

Time 

Exposure 

Time 

(sec) 

Closest D3 
Filtergram 
Time 

24 

15:46:07 

9 

15:46:00 

25 

15:53:20 

2 

15:53:30 

26 

15:54:00 

2% 

15:54:00 

34 

15:57:05 

% 

1 

15:57:00 
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Plate 2k has an exposijre time of 9 sec. and does not appear usually to he 
greatly solarized. Upon calculating the intensity at the 625OO points in the 
scanned area ve find however, that 2l6o8 are overexposed and do not fall on 
the established H & P curve. This means that only a lower limit to the 504 a 
intensity can he established over the entire flare area. Plates 25 and 26 
have about $00 data points exceeding the limits of the H & D cvirve. All of 
these data points are overexposed. The overejq)osed area is about 15" x 50" 
with the long axis running W and S. This area corresponds very closely to 
the D, flare region. Fig. II. C. 4 shows the relative intensity contours for 
URL plate 26. The average background intensity is shown on the Pig. The 
intensities have been averaged over 10.4” x 10.4" areas. The area shoxm in 
the Fig. is 155" x 135". Intensities are shown in terms of the maximum 
intensity of 100. The larger 90 contour represents the flare region. Fig. 
II. C. 5 shows this region on a scale expanded by 5 times. The entire flare 
region lies within the 100 contour region. The area inside the contour 100 
is overexposed. Figs. II. 0.6 and II. 0.7 show the 256 and 245A contour plots 
of the same region. The 256 and 245A images are not solarized. We do not 
have definitive positions of the limb for these two wavelengths and thus can 
only postulate that the intense regions here correspond to the intense region 
in Fig. II. 0.5. 

The normal exposure times for S082A were I60 sec., 40 sec., 10 sec. and 2.5 
see. Fortunately frame 34 was taken at the unusually short es^osure time of 
.31 sec. Fig. II. C. 8 shows intensity contours for this frame. Even at this 
exposure time the flare center is solarized and we cannot obtain quantitative 
intensities within the relative intensity 100 contour. 

SO82B observations were taken of the 9 Aug. 75 flare in Hell. Fig. II. C. 9 
shows the flare in the H-a line of Hell at 162OA. Unfortunately quantitative 
information was not available in time for this study. 




Figui« Il.C.if 504 a Line Intensities Relative to Maximum Intensity of 100. 

UT 15:5^J-;00 
Area 2’15'’ x 2*15" 
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Figure n.C.7 Relative 2h'^l Line Intensities 
UT 15:5^:00 
Area 29" x te" 
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S082 B POINTING 
AUG 9, 1973 IV12 FLARE 
IN AR 185 

Ha AT 15:54:00 UT 



III. A 


Theoretical Model 


The energy level model ve use consists of 19 levels in Hel^ 10 in He II 
and 1 in He III. All terras through n = 4 are included. The Tables in 
Appendix A show the energy levels and allowed lines in the model. 

A plane -parallel slab irradiated on one side by a Planckian photo spheric" " 
radiation field at radiation tenperature T^ was chosen as the geometric model. 
The gas is assumed to be optically thin for all lines and continuua xinless 
otherwise specified. 

The rate'equations for a statistically steady state have been set up as 
discussed in Rlj Kulander, 1965j Kulander, 19?6. The solution of the rate 
equations for the level populations is discussed in Appendix B. All radiative 
and collisional transitions between all terms were considered with the 
exception that excited levels only produce gro\ond state ions upon ionization. 
These equations include an external Planckian radiation field at temperature 
T^, and an adjustable parameter multiplying each line and continuum radiative 
intensity. The cross sections and reaction rates used have been discussed in 
R1 and Benson and Kulander^ 1972 • Solutions can be obtained parametrically 
for any arbitrary line and continuum radiation fields for specified election 
temperatirre and density. 

The intensity in the given line from a plane-parallel layer is determined 
by the radiative transfer equation. 
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III.A.l 



•where cos p. represents the angle between the direction of propagation and 
the outward normal = J is the source function^ is the 

linear absorption coefficient and the specific intensity of the radiation. 
The source function on a microscopic basis can be written as 

2hv^ 1 

r, / S7 — 7 -^.’ III.A.2 

where g represents the statistical weight. The populations n^ and n^ are 
obtained from solution of the statistically steady state rate equation which 
for state i is 


s 

j=i 


[n.P.. 

0 03 . 


n.P. .) = S n.P.. = 0 

X ^3 QX 


III. A. 5 


P.. 


IX 


s 

OT^i 


P. . 


where P. . is the total transition rate from i to j per second per particle 
^ J 

in the i state. 

As shown in appendix B the line source function may be written 



r + e + 

UJ0 


1 + e + T] - 6i 


III. A. 4 


where is the mean intensity^ is the Planck function at T^, is the 

normalized line profile i* = i/B and e, n and i* are non-LTE parameters 
obtained from the rate equation matrix. Tj and i are specified in Appendix B. 
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III.B.l Non -Equilibrium Parameters and Optical Depths 

To obtain estimates of the line intensities we first need to know the non- 

equilibrium parameters e , T] , and i that enter the source function. 

e is simply the ratio of collisional to radiative de-excitation of the upper 

* 

line level- 7] and i must be obtained from the full rate equation matrix. 

We have obtained these values for the 584, 304 A and D., lines for electron 

temperatures 10 through 5 X 10 ° K and electron densities 10 through 

10^^ cm~^. Tables III.B.l through III.B.5 illustrate e , T1 , and i 

for selected n and T values for the eases T, 1, 3 and 8 to be discussed 
e e 

shortly. (The number following each entry is the power of 10 by which the 
entry is multiplied.) 

* , 

For the D^ line i /T1 is shown in Fig. III.B.l along with the theoretical 
limiting value 

= 1 - 

which obtains in LTE . 

The values given here are constant with depth since the various rates entering 
the rate equation matrix have been assumed constant. This is not in general 
true even for a layer of constant electron temperature since the various radiation 
fields are functions of position. We may obtain an indication of the effects of 
these varying radiation intensities by solving for the non-equilibrium parameters 
with different assumed radiation fields. Basically these radiation fields are 
functions of optical thickness becoming dominated by the external radiation- when 
thin and approaching equilibrium values when completely thick. 

In El we made a parametric study of the effect of optical thickness in the 
resonance lines and continua upon the level populations . Level populations were 
obtained for various physically meaningful combinations of resonance lines and 
continua being optically thick or thin. Each resonance line was thus assumed to 
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TABLE III.B.l 

CHABACTEBISTIC VALXJES OP e, % i* 


\ _ « Line Case 

T^(°K) (cm~^) e T[ l'' (A) # 


10^ 10^® 

4.1-9 

2.3-4 

2.9-18 

584 

T 


It 

tf 

5.1-14 

M 

3 


1.2-8 

1.2-8 

X# 8-28 

3o4 

T 


tl 

It 

ft 

tt 

3 


4.0-8 

2.7-4 

9.4-7 

584 

T 


It ^ “ 

It 

tf 

ff 

1 


ft 

It 

9.5-3 

ft ’ 

8 


2,3o4 

6.1-2 

5.8-3 

5876 

T 


n 

5.0-2 

8.8-3 

tf 

1 


If 

5.3-2 

7.8-3 

It 

8 


1.2-7 

5.8-8 

4.8-8 

304 • 

T 


ff 

tf 

4.8-8 

11 

1 

10« 

tl 

tf 

5.2-2 

It 

8 

4.0-6 

1.5-3 

9.6-5 

584 

T 


ff 

Tf 

ft 

ft 

1 


ft 

It 

1.8-2 

It 

8 


2.3-4 

1.02 

5.7-2 

5876 ■ 

T 


If 

.84 

.10 

It 

1 


11 

.49 

.18 

If 

\ 8 


1.2-5 

4.8-6 

4.7-6 

304 



ft 

tt 

ft 

tt 

1 


It 

Jt 

4.7-2 

ff 

8 

4.0-5 

8.0-3 

4.0-14 

584 

T 


ft 

It 

8.3-13 

11 

3 


1,2-4 

4V8-5 

1.8-24 

304 . 

T 


f 1 

tt 

ft 

ft 

3 


1j6 



TABLE III.B.2 


CHARACTERISTIC VALUES OP e, T], 


t^Ck) 

""‘3 

(cm 

€ 

T) 

•55- 

i 

Line 
. (A) 

Case 

JL 

IT 

2 X 10^ 


1.1-8 

2.3-4 

6.1-13 

584 

T 



1 ! 

ft 

5.8-3 

It 

3 



8.2-9 

1.2-8 

2 . 5-18 

3o4 

T 


10^" 

If 

fl 

tr 

n 

3 


1.1-7 

2.6-4 

8.3-7 

584 

T 



tf 

It 

ft 

If 

1 



fl 

It* 

1.7-2 

ft 

8 



2.8-4 

5.4-2 

3.4-3 

5876 

T 



ff 

6.0-2 

2.8-3 

It 

1 



tl 

6 . 9-2 

1.9-3 

n 

8 



8.2-8 

5.5-8 

3.^-8 

304 

T 



11 

IT . 

11 

n 

1 


10^5 

n 

It 

6.3-4 

If 

8 


1.1-3 

2.8-3 

7.1-5 

584 

T 



n 

IT 

It 

n 

1 



ft 

If 

4.2-2 

It 

8 



2.8-2 

1.19 

4.2-2 

5876 

T 



fl 

.92 

5.9-2 

II 

1 



n 

.55 ■ 

8.1-2 

ft 

8 



8.2-6 

5.7-6 

3.^-6 

3o4 

T 



fl 

If 

1.2-5 

It 

1 


l4 

10 

11 

11 

1.5-3 

ft 

8 


1.1-4 

2.2-2 

6.4-9 

584, 

T 



It 

It 

^.9-7 

ft 

3 



8.2-5 

5*. 8-5 

2.5-14 

304 

. T 



ft 

tl 

tf 

ft 

3 


^7 





03ABIE III.B.3 


CHAEACTEKESTIC VALUES OF e, T], i* 


T Ck) 

(cm 

e 

T \ 

i 

Line 

.(A) 

Case 

JL 

TT 

3 X 10 . 

• io"° 

1.9-8 

2.5-4 

4,3-11 

584 ■ 

T 


. 

IT 

II 

IT 

It 

3 



6.7-9 

1.4-8 

6.1-15 

3o4 

T 



n 

tl 

4.2-9 

It 

, 3 


10^ 


(i.5o4) 

(9.8-7) 





1.9-7 

2.7-4 

9.2-7 

584 

T 



tt 

11 

It 

n 

1 





(4.2-2) 





it 

tr 

4.0-2 

ITT 

8 




(6.5-2) 

(1.5-3) 





3.0-4 

7.5-2 

1.5-3 • . 

5876 

T 



It 

tt 

n 

IT 

1 




(5.7-2) 

(2,2-3) 





If 

6.6-2 

2.1-3 

II 

8 




(7.0-8) 






6.7-8 

7.9-8 

3.0-8 

304 

T 



n 

IT 

11 

II 

1 



11 

It 

5.0-3 

tt 

8 


10^^ 


• (5.7-3) 

(2.2-5)' 





1.9-5 

4.1-5 

7.5-5 

384 

T 



n 

1! 

IT 

* IT 

1 



M 


(6.9-2) 





IT 

ft 

6.8-2 

n 

8 




( .88) 

. 





3.0-2 

1,30 

5.6-2 

5876 

T 



It 

tt 

• M 

IT 

1 




( .68) 

(5.6-2) 





It 

;68 

5.7-2 

tt 

8 




(1.1-5) 


\ 




6.7-6 

1.1-5 

3.0-6 

304 

a? 



U 

tt 

n 

n 

1 



tt 


7.0-5 

II 

8 


10^^ 

6.7-5 

1.2-4 

5.8-11 

584 

T 



II 

It 

1.4-8 

tl 

3 



1.9-4 

3.4-2 

4.0-7 

3o4 

a? 



II 

It 

It 

IT 

3 



,TABLE 


CHAEACTERISTIC VM.UES OF e, % i* 


Tg(°K) 

»o 

(=m-5) 

e 

Tl 

i 

Line 

(1) 

Case 

X 

7T 

5 X 10^ 

io“ 

5.8-8 

2.5-4 

1.6-9 

• 584 

T 



II 

2.0-4 

9.0-10 

It 

5 



5.2-9 

1.9-8 

3.4-12 

3o4 

T 



IT 

11 

8.1-6 

II 

5 



5.8-7 

2.8-4 

1.2-6 

584 

T 



IT 

• 2.5-4 

9.2-1 

tl 

1 



It 

M 

3.6-2 

n 

8 



3.h-k 

8.1-2 

1.2-3 

• 5876 

T 



n 

6,4-2 

1.1-3 

It 

1 



It 

4.0-2 

II 

II 

8 



5.2-8 

1.5-7 

3.0-8 

304 

T 



II 

It 

II 

II 

1 



II 

II 

1.2-2- 

It 

8 



5.8-5 

6.1-3 

8.6-3 

584 

T 



IT 

5.9-5 

II 

IT 

1 



II 

6.5-5 

.12 

IT 

8 



3.^-2 

1.43 

2.8-2 

5876 

T 



n 

l.4o 

II 

II 

1 



It 

.91' 

4.0-2 

It 

8 



5.2-6 

2.6-5 

2.7-6 

304 

T 



n 

n 

It. 

It 

1 



If 

II 

2.5-2 

It 

8 


10^^ 

5.8-ii- 

2.7-4 

3.0-8 

504^ 

T 



tl 

5.8-2 

II 

II 

5 



5.2-5 

2 .‘7-4 

3.0-8 

3o4 

T 



II 

IT 

2.7-5 

It 

5 
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have a net radiative bracket (NEB) of either 0 (completely thick) or 1 (thin). 

For each specified electron temperature, calculations were made for combinations 
of net radiative brackets corresponding to layers of varying total thickness. 
There are three resonance lines in the model for both Hel and Hell. These 
six lines together with the Lyman continue for each ion are allowed to become 
optically thick in our calculation. Thus, there are eight lines and continua 
which can be optically thick or thin depending upon the physical thickness of 
the layer. 

Nine solutions for the level populations were made for each and n^ . 

These refer to a solution where all lines and continua were optically thin 
(labeled T) and eight other solutions labeled 1-8 which correspond to the 
combinations of optical thick lines and continua given in Table III.B.6. In 
this Table (l), (2) and (3) refer to the 537 and 522 A lines of Hel 

respectively, and (U) refers to the Lsnnan continuum of Hel. Similarly the 
numbers (5-8) refer to the 30l<-, 256, 2^-3 A lines and the L 3 ?man continuum of 
Hell respectively. When a number appears in Table III.B.6 the corresponding 
transition has been assumed optically thick (i.e., in radiative detailed balance) 
in obtaining rate coefficients for the level population solutions. Transitions 
not appearing in the Table are assumed optically thin. It is noted that 
progressing from Case 1 to 8 corresponds in general to the layer becoming 
thicker. Case 8 always corresponds to all eight lines and continua becoming 
optically thick. For example, with T^ = 30,000° K Case k represents the 
NEB - 0 in the first three resonance lines of Hel and the first resonance line 
of Hell, 

In RI we also determined the optical depths in a number of lines of a 1000 km 
thick layer. Since these results are very pertinent here we shall reproduce 
them again. For each of the same cases given in Table III.B.6 we, have obtained 
the line center optical thickness in the 584, 304, 10830 A, lines and at 
the threshold of the Lyman continua of Hel and II. 
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The optical thicknesses are given in Tables III.B.7 through III. B. 11. The 
number following each entry is the power of 10 by which the entry is multiplied. 
We note that the IO83O and lines can become thick for high electron densities 

even at 50j000° K. These lines do not become thick at 10,000° K. There are many 
cases in which a number of lines and continua are optically thick. This does not 
mean, however, that simultaneous transport equations must be solved for these 
lines and continua. Which line and continuum radiation fields must be obtained 
simultaneously depends upon the level population being sought as well as the 
temperature, densiiy, and layer optical thickness. 
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Table III.B.? 

OPTICAl. OHICKWESS - 1000 KM lAYER 

T = 10,000 ®K 
e 


Case 


TRANSITION 

n 

e 

T 

I 

2 

5 

4 

5 

6 

7 8 


10^" 









584 


5.7+i^* 

= 

= 

= 

1 . 5+4 

= 

= 

= = 

504 


7.8 


= 

s 

4.1 

=: 

= 


5876 


2 . 8-8 

1 . 1-5 

2 . 1-5 

2.7-5 

8 . 0-4 

= 

= 

= = 

10830 


2.6-7 

1 . 0-4 

1 . 9-4 

2.5-4 

7.4-3 

- 

t= 


304 


1.3-2 

9.6 



8 . 6+2 

=: 

=r 


228 

10^^ 

1 . 4-6 

9.5-4 

1 . 8-3 

2 . 4-5 

8.5-2 

= 

“ 


' 584 


5.7+5 

= 

= 

= 

= 


- 


504 


76 

= 

= 

- 

= 


= 

= = 

5876 


1.3-6 

1.7-4 

2 . 2-4 

5.7-4 

6 . 7-3 

- 

=: 

= = 

10850 


1.2-5 

1.2-5 

2 . 1-5 

3.4-3 

6 . 2-2 

- 


- = 

504 


.046 




1 . 6+3 


1 

= - 

228 

10^ 

4 . 4-6 

1 . 0-3 

1.9-3 

5.1-3 

.15 


= 

= = 

584 


3.7+6 

= 

= 


= 


= 

= 

504 


7.6+2 

= 

= 

= 

= 


= 

= = 

5876 


1.5-5 

1 . 1-3 

2 . 4-3 

4.7-3 

1 . 5-2 


‘ 

= = 

10830 


1 . 2-4 

8 . 6-3 

2 . 0-2 

5.7-2 

0.12 

- 


- = 

504 


0.10 



■ 

1 . 5+3 

= 

= 

- = 

228 

0 

H 

9.5-6 

1.0-3 

2 . 8-3 

6 . 0-3 

1 . 2-1 

= 

= 

= — 

'584 

1 

5.7+7 

= 


= 

- 


- 

= - 

504 


7.6+5 

= 

= 

= 

s 


= 


5876 


2 . 1-4 

1 . 0-2 

5.1-2 

4 . 1-2 

6 . 2-2 


= 

= = 

10830 


1 . 4-3 

4.7^2 

0.13 

0.17 

0.26 


= 

- = 

304 


0.56 


1.3+22 

2 . 5+2 

1 . 1+5 

- 

= 

= rc 

228 

io"'‘ 

5.4-5 

2.7-5 

1.3-2 

2 . 4-2 

1 . 0-1 



i ^ 

584 


3.7+8 

= 

= 

= 



=: 

= 

5 o 4 


7.6+4 

= 

=: 

= 



= ■ 

- = 

00 

LT\ 


5.6-5 

0.15 

0.24 

0.26 

0.29 

= 


= = 

10830 


1 . 4-2 

0.37 

0.54 

0.59 

0466 

- 


- = 

304 


4.7 

1.7+2 

5.0+2 

6 , 1+2 

9.5+2 

— 


' = = 

228 4 . 5-4 1.6-2 4 . 8-2 5.9-2 

* The nmnber folloi-jing each entry is the power of 

9.0-2 

10 by i^hich the entry 

is multiplied.. 
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Table 


OPTICAL THICifflESS - 1000 KM lAIER 

T = 2 X 10^ ° K • 
e 


Case 


TRABSITIOW 

n 

e 

T 

1 

2 

5 

4 

5 

6 

7 

8 


10 ^*^ 










584 


1.7+5 

5.4 

0.92 

0.52 

= 

- 

= 

7.0-5 

= 

5 o 4 


0.36 

7.2-4 

1 . 9-4 

1 . 1-4 

- 

= 

- 

1.5-6 


5876 


1.7-5 

1 . 8-5 

1.7-3 




= 

9.6-4 

= 

10850 


1.6-2 

1.7-2 

1.6-2 


= 



9.5-5 

= 

304 


1.5+5 

2 . 8+5 

= 

= 

= 


= 

2 . 0+3 

= 

228 

10"" 

0.10 

0.19 




= 


0.13 


584 


2.0+4 

2 . 2+2 

61. 

26. 

= 


= 

1.0 

= 

504 


5.7 

6.3-2 

1.7-2 

7.4-3 

= 

= 

r: 

2 . 9-4 

= 

5876 


4.5-2 

0.11 

= 

= 

= 

= 

- 

5.5-2 

= 

10850 


0.38 

1.0 

0.86 

0.95 

= 



2.9-4 

— 

5 o 4 


5.1+5 

2 . 1+4 

= 

= 


=5 

= 

= 

= 

228 

10 ^ 

0.45 

1.8 

- 


= 

= 


2= 

= 

584 


2.3+5 

1 . 0+4 

2.0+3 

5.6+2 

= 


= 

100. 


504 


65. 

2.8 

0.56 

0.16 

= 

= 


2 . 8-2 

= 

5876 


0.90 

4.3 

4.0 

3.1 

= 


= 

1.0 

= 

10830 


5.3 

24 . 

22 . 

17. 

= 

= 

- 

= 

= 

5 o 4 


5.3+4 

2 . 0+5 

= 

= 



c:: 


= 

228 

H 

0 

2.7 

16. 

=: 

- 



= 

= 

= 

584 


2 . Oi '6 

1 . 2+5 

1.7+4 

8.5+5 

= 

= 

= 

5.2+3 

= 

5 o 4 


6.1+2 

57. 

5.2 

2.6 

= 



1.6 


5876 


18. 

65. 

48 . 

58. 

= 


= 

28 . 


10830 


57.. 

1.3+2 

90.. 

69. 

= 

= 

= 

52. 


304 


4.0+5 

1.9+6 



= 


= 

= 

= 

228 

l 4 

10 

56. 

170. 

~ 


= 

= 


1 


584 


2.0+7 

1 . 2+6 

3 . 8+5 

3.0+5 



- 

2.6+5 

- 

504 


5.8+3 350, 

110. 

87. 

- 



75. 


5876 


2 . 8+2 

7.7+2 

5.4+2 

5.1+2 



- 

4 . 8+2 

= 

10830 


3.5+2 

8 . 8+2 

5.&+2 

5.4+2 

= 



5.1+2 


304 


5.1+6 

1.9+7 

= 

=: 


- 

- 

— 

- 

228 


450. 

1.7+5 

=: 

56 “ 




- 

= 


56 



Table III.B.9 


OPTICAL TEECKKESS - 1000 KM LAYER 

T = 3 X 10^ ° K 
e 


Case 


TRANSITION 

n 

e 

T 

1 

2 

5 

4 

5 

6 

7 

8 


10"° 










58U 


23. 

= 

22. 

21. 

2.5-2 

5.2-3 

2.7-3 

1.7-5 

1.6-5 

5 o 4 


8.4-3 


- 

8.2-3 

1.0-6 

2,0-7 

1.0-7 

5.9-6 

6.0-11 

5876 


1.5-3 


= 

1.2-3 

7.7-4 

7.0-4 

= 

4.3-6 

3.2-6 

10830 


1.2-2 

= 

= 

1.1-2 

7.2-5 

6.5-3 

6.4-3 

4.0-5 

5.0-5 

304 


1.6+3 

= 

= 



= 

= 

10. 

= 

228 

10"" 

0.18’ 



0.17 

= 

r:: 

= 

1.1-5 


584 


6.3+2 

= 

= 

= 

2.0 

0.43 

0.15 

8.8-5 

2.6-4 

504 


0.19 

= 

= 


7.2-4 

1.5-4 

5.0-5 

5.2-6 

1.0-7 

3876 


6.8-2 

= 


“ 

5.2-2 

4.5-2 


2.9-3 

8.2-5 

10830 


0.57 

= 



0.43 

0.37 

= 

2.4-2 

7.0-5 

304 


1.6+4 

= 

= 

= 

= 

- 


1.0+3 

= 

228 

10^ 

1.7 

r 


= 


= 

== 

0.10 


584 


7.5+5 

- 


- 

62. 

9.6 

2.5 

1.0 

0.13 

504 


2.6 



=: 

.22 

5.4-3 

8.9-4 

3.6-4 

4.6-5 

5876 


1.5 


- 

=: 

1.4 

= 

1.1 

0.42 

0.15 

10830 


6.8 

= 

= 

= 

6.2 

5.9 

4.7 

1.8 

0.65 

304 


1.6+5 


= 

= 

= 

= 

= 

6.4+4 


228 

10« 

18. 





— 


7.0 


584 


4.8+4 

= 

= 

= 

6.5+2 

76. 

38. 

32. 

17. 

504 


17. 




0.23 

2.7-2 

1.4-2 

1.2-2 

6.4-5 

5876 


22. 

“ 


=: 

18. 

15. 

12. 

10. 

7.5 

10830 


32. 

= 

= 


25. 

19. 

15. 

15. 

10. 

304 


1.6+6 

= 

= 


= 


— t 

1.4+6 

- 

228 


1.8+2 

= 



= 



1.5+2 

= 

584 


5.1+5 

5-3+5 

3.7+5 


5.9+5 

1.8+5 

1.4+5 

= 

1.1+5 

■504 


1.2+2 

= 

1.3+2 


2.5 

0.69 

0.48 

= 

0,38 

- 5876 


2.0+2 

- 

2.4+2 


1.9+2 

1.4+2 

1.3+2 

= 

1.2+2 

10830 


2.6+2 

= 

2.1+2 

= 

1.7+2 

1.2+2 

1.1+2 

- 

1.0+2 

304 


1.5+7 

= 

=5 

= 

1.6+7 

= 

= 

= 

— 



Table III, B. 10 


OPTICAL THECKIffiSS - 1000 KM LAYER 

T = 4 X 10^ ® K 
e 


Case 


TRANSITION 


T 

1 

2 

5 

4 

5 

6 

7 

8 


.10 

10 










584 


1.9 

= 

1.5 

0.37 

1 . 9-4 

1 , 0-7 

1.7-8 

8.5-9 

5,5-11 

504 


7.6-4 

= 

4.5-4 

1.5-4 

7.6-8 

4 . 0-11 

6.8-12 

3 . 4-12 

2 . 2-14 

5876 


8.5-4 


6 . 4-4 

1 . 6-4 

1 . 6-8 

1 . 4-8 

1 . 2-8 

= 

6 . 6-9 

10850 


7.7-3 


5.9-3 

1.5-3 

7.7-7 

1.3-7 

1.0-7 


6 . 0-8 

304 


1.5+3 

= 

1.2+5 

2.9+2 

2.9-2 


= 


= 

228 

10^^ 

0.18 


0.11 

.035 

3.5-6 

= 


— 


584 


4 o. 

58 . 

33. 

15. 

6.7-3 

4.5-5 

8 . 0-6 

2.7-7 

4.7-9 

5 o 4 


1.7-2 

= 

1.5-2 

7.0-5 

2 . 8-6 

1 . 8-8 

5 . 4-9 

1 . 1-10 

2 . 0-12 

5876 


0.52 

0.51 

0.27' 

0.15 

5.6-6 

= 

5.8-6 

- 

1 . 2-6 

10850 


0.26 

0.25 

0.21 

0,11 

4.5-5 

- 

5.9-5 

= 

1.5-5 

504 


1 . 4+4 

= 

1.3+4 

7.6+3 

2.8 

= 


= 

= 

228 

10^ 

1.8 


1.6 

0.81 

3 . 6-4 

— 

= 



584 


4 . 2+2 

k -4 1+2 

3.7+2 

2,5+2 

0.82 

8.1-3 

1.1-3 

2.5-4 

3.3-5 

5 o 4 


0.17 

- 

0.15 

0.11 

5 . 4-4 

5.4-6 

4.7-7 

1 . 2-7 

1 . 4-8 

5876 


0.69 

■ 0.67 

0.61 

o. 4 i 

1 . 4-5 

= 

- 

1 . 1-3 

4 . 0-4 

10850 


2.6 

2.5 

2.3 - 

1.5 

5.0-3 

= 

4.9-5 

4.0-5 

1 . 4-5 

504 


1.5+5 

1 . 4+5 

1 . 3+5 

8.7+4 

2.9+2 

= 

= 

= 

= 

228 


18. 

17. 

16. 

11. 

3 . 6-2 

= 

= 

= 



loW 



; 







■ 584 


2.7+3 

2.5+3 

2:3+3 

1.6+5 


0.80 • 

0.10 

4.5-2 

2 . 0-2 

504 


1.0 

0.92 

0.84 

0.68 

2 . 2-2 

2 . 1-4 

3.6-5 

1,8-5 

7.7-6 

5876 


9.1 

8.7 

8.0 

5.6 

0.19 

0.16 

0.13 

0.10 

8 . 0-2 

10850 


11 . 

10 , 

9.2 

6.5 

0.22 

0.19 

o.i 4 

0.12 

8.5-2 

504 


1,5+6 

1 . 4+6 

1 . 2+6 

8 . 8+5 

3.0+4 



= 


228 


180. 

170. 

160, 

100 , 

3.6 

= 



= 

584 


2 . 1+4 

2 . 0+4 

1.7+4 

1 . 2+4 

3.1+3 




9.0 

504 


8.4 

7.8 

6.4 

3.0 

0.77 

1.5-2 

4 . 0-5 

3 . 2-3 

2 . 2-3 

'5876 


99.. 

95. 

83. 

58. 

15. 

12. 

10, 

9.0 

8.5 

10850 


76.* 

72; 

62.' 

45: 

11 . 

10' 

7.0 

6.6 . 

6.0 

504 


1 . 4+7 

= 

1 . 2+7 

9.0+6 

2.3+6 

~ 


- 

- 

228 


1.8+5 

1 - 7+3 

1.4+5 


5.0+2 

= 

= 

= 

= ' 



Table III. B. 11 

OPHCAL THICKHESS - 1000 KM lAYHK 

h 

T = 5 X 10 °K 
6 


Case 


TRANSITION 


T 

1 

2 

3 

4 

5 

6 

7 

8 


0 

1 — 1 

0 

H 










584 


0.55 

0.21 

4 . 2-2 

5.5-5 

2.1-7 

- 




504 


1 . 5-4 

9.0-5 

1,8-5 

1.5-6 

9.0-11 

- 

- 


- 

5876 


5.4-4 

3 . 4-4 

6.5-5 

5.3-6 

5.2-10 


- 

- 

~ 

10850 


4.9-3 

5.0-5 

6 . 0-4 

5.0-5 

3.0-9 

- 

- 

- 

- 

504 


1.5+3 

7.8+2 

1 . 6+2 

15. 

7 . 8-4 

= 

= 

- 

= 

228 

10^ 

0.18 

0.11 

2 . 2-2 

1 . 8-3 

1 . 1-7 

— 




584 


6.6 

4.0 

1.5 

0.26 

4 . 0-5 

- 

- 

- 

- 

5 o 4 


5.0-3 

1.8-5 

5.9-4 

1 . 2-4 

1.8-8 

- 


- 

- 

5876 


1.8-2 

1.1-2 

5.6-5 

7.2-4 

1 . 1-7 

- 

- 



10850 


o.i 4 

8 . 4-2 

2 . 8-2 

5.6-5 

8 . 4-7 

-- 

- 

- 

- 

504 


1.3+4 

7 - 8+3 

2 . 6+5 

5.2+2 

7.8-2 

= 


“ 

- 

228 

10^ 

1.8 

1.1 

0.36 

7.2-2 

1 . 1-5 


— 

“ 

“ 

584 


62. 

37. 

12 . 

5.1 

5.7-3 

- 


- 

- 

504 


5.0-2 

1 . 8-2 

5.8-5 

1.5-5 

1 . 8-6 

- 

- 

- 


5876 


0.37 

0.22 

7,4-2 

1 . 8-2 

2.2-5 

- 

- 

- 

- 

10830 


1.2 

0.75 

0.24 

6 . 0-2 

7.2-5 

- 

- 

- 

- 

504 


1.5+5 

7 . 8+-4 

2 . 6+4 

6.5+3 

7.8 

= 


== 


228 

loW 

18. 

11 . 

3.6 

0.90 

1 . 1-5 





584 - 


4 . 1+2 

2 . 4+2 

82. 

20 . 

0.25 

4 . 0-5 

4.5-4 

2 . 3-4 

9.0-5 

504 


0.18 

0.11 

5.6-2 

. 8 . 8-3 

1 . 1-4 

1 . 8-6 

2.0-7 

1 . 0-7 

4 . 0-8 

5876 


4.8 

2.9 

0.96 

0.24 

2.9-5 

2 . 6-5 

2 . 2-3 

1.7-5 

1.3-3 

10850 


. 4.9 

2.9 

1.0 

0.25 

2.9-5 

2 . 6-5 

2 . 1-3 

1.7-5 

1.5-5 

5 o 4 


1.3+6 

7 - 8+5 

2 -. 6+5 

6.3+4 

7.8+2 

= 

= 

= 

= 

228 


180. 

no. 

36. 

8.8 

0.11 

— 

1 



584 


5.0+3 

1.9+3 

5 . 0 f 2 

H 

0 

19. 

0.35 

0.11 

8.5-2 

5 - 9-2 

5 o 4 


1.4 

0.84 

0.22 

J.0-2 

8.9-3 

1 . 6-4 

5.1-5 

4 , 0-5 

2 * 8-5 

15876 


47. 

30. 

7.8 

2.4 

0,29 

0.23 

0.20 

= 


10850 


33. 

21 . 

5.6 

1.7 

0.21 

0.18 

0.14 

- 

- 

504 


1.3+7 

7.8+6 

2 . 1+6 

6.5+5 

7,8+4 

= 

= 

- 

- 

228 


1 . 8^5 

1 . 1+3 

2.9+2 

9.0+1 

11 . 

= 


- 
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III.B.2 D, Line Profiles and Transmitted Intensities 

Using the line source function of Eq. III.aA for a finite plane - 
parallel layer above the photosphere we may obtain the emergent intensity 
tl (0 ^|j,)/b ] =s I ^(0,[a) from the surface of a layer of optical thickness t_ 

V 6 V i- 

at angle cos p, •vri.th respect to the normal, from 

(0,h) = J exp(- J cp^ dt/|x) dt/h * III.B.l 

o o 


With cp a constant ifith depth, the photo^heric radiation being given by B 

^ X* 


and assuming no continuous absorption the emergent intensity becomes 
(see Appendix B for details) 


I (0,p) = CO cp 2 L 


<p + k 


of* 


*1 ■ 

* * (^v-V> ■ 


+ 00 cp 2 

a a 


-k t, cp t./p 

a 1 “v p 

e - e 


-Pv - “ 


+ 4 


(l - e' " 





- ‘I’v Vt* 


III.B.2 


This is a homogeneous slab solution assuming an Eddington approximation 
which we shall use in obtaining transmitted intensities. 


B5 line profiles have been calculated for a wide range of values of 
e, T|, i and from Eg,* Figure illustrates results for 

r n p ^ 

T = 6000'^K, T = 20,000'^K, n = cm ^ for various values of e, T|, i 
re e 

and tj^. 

The filter has a Gaussian profile with a total half width of O.ho k. 
The transmission of the filter can he vjritten as 


T = 


v?here v 

% 


Exp - 

and 

^ /Sc T^/m 


/vo2. . 
0.20 ^ 


III.B-.3 
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Some characteristic values for AXj^ are given in Table III.B.l. The total 
integrated, intensity observed through the filter is 

J \ • 


¥e define the total transmission TT of the filter as 


TT = 


J \ ly 
J 


Fxgure III.B.2 illustrates values of TT for the case e = 10 Te = 20,000 K, u = 1, 
for various values of Tj and i . The curves are labeled mth.the exponents of 11 

* -5 «• _li , 

and i, for exaiiip3ethe curve T] « 10 ^ i =10 is labeled 3,4. 


The most useful quantity for our interpretation is the total transmitted 
line intensity ivith respect to the total transmitted continuum intensity, i.e. 

f T I dv 

R = Ljlj^ . 

Selected values for R are shoTO in Figures IH.B.3 to III.B.6 for T = 6O0 CPk, 

4 1 4 * -1 -2 -3 ^ -4 

|i = 1 , Te = 2 X 10 and 4 x 10 and for e values 10 ,10 , 10 and 10 

-2 -3 * - 1 }. 

The curve for e = 10 , 1] = 10 i =10 for example is labeled 2,3,4. 

has been set equal to B(T^). 
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<?' f 0.2 0.3 o. <f o. ^ 0.6 0.7 

T 

Figure UI,B.5 Total Transmission for T = 20,000° K, e = lo”^, H = 1. ^ 

Curves are labeled mth negative exponent values of T| and i . 
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Figure III.B.7 Total Transmitted Line Intensity mth Respect to Continuum for e s= 10"^. 
Curves are labeled •with negative exponent values of T] and i*. 
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III. C. Electron Temperatures and Densities 

The maximum flare D^ intensities are about 1.2 times the normal background 
continuum. For each electron temperature and density we ask whether the 
corresponding values of T| will give line profiles yielding trans- 

mitted filter intensities greater than the continuimi. In order to obtain an 

absolute D_ intensity greater than the continuum with a layer no thicker than 
^ 11 

several thousand kms^ n must be at least lO . The following temperatures, 

® 15 

densities and optical thicknesses provide sufficient intensity; n = 10 , 

0,0 o ® 5 

T « 20000 K - toOOO K, t, = 0 - 2j n « 10 , T = 20000 K, t^ = 10 - 10 . 

-L e 1 

6 2 » 

The maximum D_ flux in the flare is approximately 5*^ x 10 ergs/cm /sec/A 

^ "5 2 o 

at the sun. The transmitted intensity is 6.0 x 10^ ergs/cm /sec. At 20000 K 

the total D^ line flux is .157 times the intensity or 5*5 x 10^. 

We cannot determine the maximum 50kA intensities due to the solarized 
condition of plates 24, 25, 26 and 54. On plate 26 the 100 contour is at 40 
times the background continuum while on plate 3^ it is about Ito times the 
background continuum. By examining various positions within the flare we can 
determine that the intensities in the flare at 15:57:05 have decreased about 75^ 
from those at 15:54:00. 

The 256A line intensities on plate 26 have not solarized. If we assume 
the same ratio of increase from positions near the intensity 100 contour on 
the 304 contour map as on the 256 contour map we can estimate the 304 a central 
intensities. Ihis yields a maximum 504 a flare intensity of 80O times back- 
ground intensity. 

2 

We shall use a quiet sun background intensity in of 8600 ergs/cm^/sec/st 

T 2 ' o 

or flux 2.2 X 10* ergs/cm /sed. We thus arrive at an approximate 504 a to 

flux ratio of 40. Using the intensity curves derived in Rl we can determine 

the temperatTjres and densities allowing for a 504 a/D_ ratio between 20 and 60. 

4^0 12 15 -5 

We find the following acceptable ranges: T = 3 x 10 K, n =10 - 10 cm 

, 4 o 12 12 “5 ® 

T = 4 X 10 K, n = 10 - 3 X 10 cm . In order to have the correct 

® 4 

absolute intensities the following optical thicknesses are required: T = 3 x 10 , 

12 . , _ _ „ „4 _13 . , _4 


n =10 , t = .1 - .3,- T = 3 X 10 , n = 10 t, = .1 - .2; T = 4 x 10 , 

® 12 ^ 46^1 

n =10 , t = .1 - .5; T = 4 X 10 , n = 3 x 10 


g , ^ ' e ' , t^ = .2 - . 5 . Comparing 

these ranges with those consistent with the D, emission we find the most 

5 4 o 

probable flare condition at the brightest point to be T ~ 3 x 10 K, n 


10« om-5 


and t^ 


.1 - . 2 . 


Further inte3:pretation is undervra-y. 
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Appendix A 
EKERGY LEVEL MODEL 


TABLE I 

He Ij II Energy Levels 


He I - i = 1 


Energ^f (ev) Have Nos. 


1 

T 2 

Is 

"s 

0 

0 

2 

ls2s 


19.821 

159850 

3 

ls2s 


20.618 

166272 

4 

ls2p 


20.966 

169081 

5 

ls2p 


21.220 

171129 

6 

ls 3 s 

h 

22.721 

183231 

T 

ls 3 s 

h 

22.923 

184859 

8 

ls 3 p 

'1 

23.009 

185559 

9 

ls 3 d 

3 d 

23.076 

186096 

10 

ls 3 d 

h ) 

23.076 

186099 

11 

ls 3 p 

IpO 

23.089 

186204 

12 

Isii-s 

3 s 

23.596 

190292 

13 

IsM-S 


23.676 

190935 

m 

Is^tp 

3 pO 

23,710 

191211 

15 

Isi^d 

3 d 

23.738 

191439 

i6 

ls 4 -d 

Id 

23.739 

191441 

17 

ls 4 -f 

3 pO 

23.739 

191447 

18 

ls 4 -f 

IpO 

23.739 

191447 

19 

ls 4 p 

IpO 

23.744 

191487 

He II 

- i 

= 2 



1 

Is 

2s 

0 

0 

2 

2s 

"s 

40.8099 

329179.57 

3 

2p 

2pO 

40.8091 

329182.02 

iv 

3 s 

2s 

48.3662 

390140.76 

5 

3 p 

2p0 

48.3664 

390141.49 

6 

3 d 

^D 

48. 3665 

390142.64 

T 

4 s 

2s 

51.0113 

411476.98 

8 

4p 

2p0 

51.0114 

411477.28 

9 

4 d 

2d 

51,0115 

411477.77 

10 

4f 


51.0117 

411477.95 
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TAELE IT 


He I, II Lines 


He I 


Upper 

Level 

Lower 

Level 

Notation 

\( k ) 

A 

(l 0 ®/see) 

f 

4 

2 

1 

10830 

.1022 

.5391 

5 

1 

2 UV 

584,4 

17.99 

.2762 

5 . 

3 


20582 

.01976 

.3764 

6 

4 

10 

TO65 

.278 

.0693 

7 

5 

45 

T 281 

.181 

.0480 

8 

2 

2 

3889 

09478 

.06446 

8 

6 


4 . 30+4 

.0108 

.896 

9 

4 

11 

5876 

.706 

.609 

9 

8 


1.86+5 

1 . 28-4 

.111 

10 

5 

46 

6678 

' .638 

.711 

11 

1 

3 UV 

■537.1 

5.66 

.0734 

11 

3 

4 

5016 

.1338 

.1514 

11 

T 


7.43+4 

.00253 

.629 

' 12' 

4” 

12 

'4713 

'.106 

.0118 

12 

8 


21120 

,0652 

LT\ 

• 

13 

5 

4 T 

5049 

.0655 

.00834 

13 

11 


21132 

.0459 

.103 

14 

2 

3 

3188 

.0505 

.0231 

14 

6 


12538 

.00608 

.0429 

14 

9 


195^3 

.00597 

.0205 

14 

12 


1.09+5 

.0505 

1.21 

15 

4 

14 

4472 

.251 

.125 

\ 

15 

8 


17002 

.0668 

.482 

15 ' 

14 


4.39+5 

4.16-5 

.200 

16 

5 

48 

4922 

.202 

.122 

16 

n 


19089 

.0711 

.647 

IT 

9 


18688 

.139 

1.02 

IT 

15 


1 . 43+7 

6 . 01-10 

.0033 

18 

10 


18699 

.138 

1.01 

18 

16 


1.67+7 

4 . 34-10 

.00253 

19 

1 

4 UV 

522.2 

2.46 

.030 

19 

3 

5 

3965 

.0717 

.0507 

19 

T 


15088 

.0137 

.140 
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TABLE II (Continued). 



He II 


3 

1 

303.80 

100 . 

.4162 

4 

3 

1640.5 

1.01 

.01359 

5 

1 

256.3 

26.8 

,.07910 

5 

2 

1640.4 

3.59 

.4349 

6 ’ 

3 

1640.4 

10.3 

.6958 

7 

3 

1215.1 

.413 

3 . 045-3 

7 

5 

4687.0 

.294 

.03225 

8 

1 

243.03 

10.9 

.02899 

8 

2 

1215.1 

1.55 

.1028 

8 

4 

4686.8 

.491 

.4847 

8 

6 

4687.2 

.056 

.01099 

9 

3 

1215.1 

3.30 

.1218 

9 

5 

4686.9 

1.13 

.6183 

10 

6 

4687.1 

2.21 

1.018 
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Appendix B 
BASIC EQUATIONS 


In the statistically steady state the rate equation describing the 
population n^ of the state i is 

S (n.P.. - n.P. .) = s n.P.. = 0 


P. . = - S P. o 


A,1 


where P. . is the total transition rate from i to j per second per particle 

X J 

in the i state. In general, P . . = R. . + C. ■where R. . and C. . represent 

ij ij xj 

the radiative and collisional transition rates respectively. We shall 
assume a Maxwellian distribution for the electrons and helium particles and 
since 'we also assume a kno'wn external radiation field, the transitions 
involving the continuum can be represented by a single term in Eq, A.l. We 
can characterize the system of linear equations A.l by a matrix whose co- ^ 
efficients a. . are equal to P... In representing matrix elements and. 
co-factors thereof, ve shall always let the first subscript refer to the 
row and the second to the column. Because of the definition of the transi- 
tion rate the subscripts of the P* s will be reversed when substituted for 
the matrix elements a. The general solution of Eq, A.l is 


n. = X P^ j 


ra 


W 


X = 
m 2 


A.2 
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■where is the co-factor of the element a . = P. in the coefficient ma- 

Jno 

trix represented by the i equations of type A.l and IT is the total number 

3 

of helium particles per cm . The matrix of coefficients has the property 
that the co-factor of all the elements in a coliimn are equal, i.e. is 
independent of m. 


TThen the medium becomes optically non-thin for certain frequencies 
the radiation field producing internal excitation and ionization for these 
frequencies is no longer merely the external radiation field but is partly 
dependent on the internal properties of the gas and must be determined from 
the radiative-transfer equation 


V- 



A. 5 


■vdiere cos"^M. represents the angle between the direction of propagation and 
the outifard normal z and t = fk dz, S is the source function, k is 
the linear absorption coefficient and is the specific intensity of the 
radiation. In LTE B = B , however in the non-LTE case S must be specified 
in terms of microscopic processes. In terms of such processes the trans- 
fer equation governing the spectral line between upper level u and lower 
level I may be ^rritten as 


- hjtijL ® hv - n B ^ilr hv + hjtk ] I 

dz '• 1 lu V u ul^v c V 

- n A , j hv - 4rtk S (T ), A.4 

u c c e ’ 

where k^ represents the continuum absorption coefficient at frequency 

and S is the continuum source function. B_ , B _ and A are the Einstein 
c iu ul ui 

transition probabilities for absorption, stimulated emission and spontaneous 
emission, and 9^ are the normalized emission, stimulated emission 

and absorption coefficients ■within the line defined such that 

OD CO OO 

J 5/v = I Ir^dv . r ]■ J^dv a. = 1. 
o o 4jt o 





The contimiAm absorption coefficient is general3^ very small compared i-rith 
the line absorption coefficient near the line center and will be neglected 
in determining the source function id-thin the line. Using the standard 
relations between the Einstein coefficients and assuming 
the source function becomes 


« _ ^ 1 . 

c^. r(g^/i3;)(n^/n^)]-l ^ 

vhere g represents the statistical weight. The minus one term in the 
denominator represents stimulated emissions. 


In evaluating the radiative excitation rate R, . for transitions 

^ 0 

letT^een bound levels the line radiation field enters as 

CO 

o 


where 

j (t) = 7;i- r I (t,h) dti) 

is the mean intensity and du) represents the solid angle. It is thus con- 
venient to formulate the transfer eq^uation in terms of rather than 
It is now convenient to separate those components involving the unkno™ 
radiation fields from the co-factors. This is done by expanding the 

determinant in terms of its co-factors Thus 


= S + .2 . 


k/l 


'lu" 


kj^lj^u 


Xk'’ 


k/a >* k/uA 


" A,6 


Actually^ may appear in many of the co-factors since the line u-1 may 
fall in the ionization continuum of some other transition. The influence 
of as wen as the line radiation in general on the bound-free radiative 
rate's ivlll be neglected. Using the standard relationship between the 
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Einstein coefficients, equation A, 6, and remembering that 

and J* J^s^dv the source fixnction may he written as 


®ul ^ul pUl 


plu + eB + t 


1 + e + 11 


A. 7 


where B is the Planck function 


B„ - P 


lu 


V e' f-^ul C 


ul 


e = 


"ul 


A 


ul 


ul “ 2 

c 


2hv ^ 
g. 


u 


, ^ p __i X) 

V ^ul , rtUl , / /t 

A^Q VuT^l 


A. 8 
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The terms entering the numerator of Eq. A, 7 each represent a- method of 
populating the upper level from the lower level. The first term repre- 
sents direct radiative excitation, the second direct collisionaJL excitation 
and the third any combination of radiative or colli sional processes involv- 
ing one or more intermediate levels in golxjg from the lower to the upper 
level. The denominator, on the other hand, consists of terms indicating 
transition paths from the upper to the lower level. All the terms are 
normalized with respect to The first term represents direct radiative 

de-excitation, the second direct collis tonal de -excitation and the third 
any indirect process going from the upper to the lower state. 
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The problem remains of finding a to use in evaluating the source 

function. One method we shall adopt is the general procedure of Thomas and 
Zirker (1961) in ea5)loying the Eddington transfer equation and replacing the 
integral in equation (A. 2) by a standard quadrature. Specifically, an n-point 
Gaussian quadrature, as discussed by Avrett and Hummer (1965), is eniployed 
over the central line region v^ere significant transfer of radiation occurs. 

As they point out, the Gaussian quadrature proves much more accurate than the 
Hermits -Gauss quadrature for finite layers. 


The transfer equation for each quadrature frequency v^, ass-uming 
. negligible continuous opacity, is 
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■where 




and Kq are the line and line center absorption coefficients and are the 
wei^t coefficients of the Gauss quadrature. 


The solution of equation (A.9) regarding e, T) and i as constants for 
a layer of optical thickness t^, is 
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where 


A = 


e + i 
e + T| - 6i 
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and 


cDs= [(l + e+ Tl“6i)] 
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The k are given by 

VpiO 

n coa. 
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JSpplication of the bo\mdary condition, which requires no incident radiation on 
the top surface of the layer, and an intensity on the bottom surface (W is 
a dilution factor) yields the following expression for the integration constants 
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The line source function (equation (A. 7) then becomes. 
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The emergent intensity [I^(0, }i)/b] = (O, p) from the surface of the 

layer at angle cos with respect to the normal, is 
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With given by equation (A. 15) and cp^ a constant -vriLth depth, the intensity 


becomes 
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This equation is the basic equation from which the total integrated line 
radiation mil be obtained. 

When the optical thickness exceeds unity, a self -reversed or apparent 
absorption region develops near the line center resulting from the rapid 
increase in the source function with depth in the layer. 
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